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ERRATA 
Eq.  (2.9)  change  ...+d(?)  sJJ,  ^^^4).,.  to  read  .  ..+d(c)sJJ'/^^ 
Add  (j,k  =  1,2,3, U)  after  last  eqs. 
Eq.  (2.2U)  change  »  a^  »  to  read  •«  a^^,  n  , 
Eq.  (2.26)  insert  "  -  "  before  the  braces. 
Eq.  (2.27)  change  »  a^  "  to  read  "  a^^  '•  . 


Eq.  (2.28)  change  "  p^  "  to  read  "  p^',  «  , 
Eq.  (2.29)  change  "  5^  •»  to  read  ••  djj',  "  . 

Page  11   Eq.  (2.36)  add  "■?;-"  as  superscript  to  square  brackets. 


Last  line  should  read  "  where  a  =  a^.  #  B  =  3^. .  5=5^.. 

line  lU  change  "  |arg(-iY4)  >  J  "  *'°  read"  |arg(-iY4)  <  ^  "  » 

In  footnote  change  "  Appendix  A"  to  read  "  Appendix  I"  . 

Eq.  (U.l)  in  second  integrand  change  "  Q"'"(n)  "  to  read  "Q^'"(^)" 

line  6  change  "  6^"^  =  0"  to  read  "  S^-^O"  . 


Eq.  (U.l6)  change 


m  I 


t 


to  read 


m 


/(4^-i)(i-,n^) 


7(4^-1)  (i-n^) 


line  2  from  bottom,  change  "  p!^-"  to  read  "  Pp*^*  " 

line  3  change  -Y^^l   to  read  -Y^^,  , 
n+1         n+1 


Eq.  (1.17)  Chang,  "a""'  «  to  read  "  a"*   '» 


n,p  n,p 

Second  line  after  Eq.  (1.18)  change  "  Eq.  (3.7)  or  (3.11)"  to 

read  Eq.  (3.1)  and  (3.3). 
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Page  32   Eq.  (1.23)  change  <'Y^'^^(z)"  to  read  ■  ¥^^^  (z)"  . 
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Abstract 


Starting  from  a  general  integral  theorem  for  spheroidal  functions, 
we  obtain  several  integral  relations  involving  products  of  spheroidal  func- 
tions by  choosing  the  kernel  in  the  integrcil  to  be  of  the  form  e"  ^  '^  Lu, 
where  u  is  a  solution  of  the  time-harmonic  wave  equation  in  spheroidal  co- 
ordinates and  L  is  an  operator  commuting  with  the  Laplacian.  Operators 
possessing  this  property  are  certain  linear  combinations  of  the  linear  ajid 
the  ang:ular  momentum  operators  and  powers  and  products  of  these  combinations. 
Several  new  integral  relations  between  associated  Legendre  fvinctions  have 
been  obtained  as  limiting  cases  of  the  integral  relations  involving  products 
of  spheroidal  functions. 
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1,  Introduction 

In  this  paper  we  obtain  a  number  of  integral  relations  involving 
products  of  spheroidal  functions  based  on  a  general  integral  relation  due 
to  Ince"- J,  We  choose  the  kernel  in  the  integral  to  be  of  the  form 
e"  "   Lu,  where  u  is  a  solution  of  the  time-harmonic  wave  equation  in 
spheroidal  coordinates  and  L  is  an  operator  commuting  with  the  Laplacian. 
Operators  possessing  this  property  are  certain  linear  combinations  of  the 
linear  and  the  angular  momentum  operators  and  powers  and  products  of  the 
combinations. 

Spheroidal  functions  have  played  an  increasingly  important  role  in 
problems  of  applied  mathematics  and  quantum  mechanics '-  -■ .  For  example,  the 
use  of  these  functions  has  made  it  possible  to  find  exact  solutions  of  the 
problems  of  diffraction  of  acoustic  and  of  electromagnetic  waves  by  circular 
discs  and  apertures'-^.  Although  a  systematic  treatment  of  the  theory  and 
applications  of  spheroidal  functions  has  recently  been  published  by  Meixner 
and  Schaf ke  ^  -^ ,  the  theory  of  spheroidal  functions  in  general  has  not  been 
as  extensively  developed  as  the  theory  of  other  special  functions  of  mathe- 
matical physics  (Legendre,  Bessel,  Laguerre,  Hermite  functions,  etc.).  For 
the  latter  functions,  whose  differential  equations  are  of  the  hypergeometric 
type,  there  exist  well-known  recurrence  relations  connecting  three  contiguous 
functions.  These  formulas  provide  a  method  of  constructing  numerical  tables. 


and  fvirnish  valuable  insight  into  the  nature  of  the  functions  themselves. 

rec 
[6] 


Whittaker '^-'  has  given  a  method  by  means  of  which  he  derived  'reciurr- 


ence  relations'  for  Mathieu  functions.  This  method  was  used  by  Sharma' 
for  Lame"  functions  and  recently  by  Marx'' -I  for  spheroidal  functions.  It 
should  be  noted  that  Whittaker's  method  does  not  yield  explicit  expressions 
for  the  coefficients  in  these  'recurrence  relations'.  Subsequently  Marx'--' 
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was  able  to  obtain  explicit  expressions  for  the  coefficients  of  the  'recurrence 
relations'  for  the  spheroidal  functions  by  making  use  of  the  expansion  theorems 
and  the  orthogonality  properties  of  these  functions. 

However,  Meixner'-  -■  has  shown  that  the  coefficients  in  the  'recurrence 
relations'  for  Mathieu  functions  given  by  Whittaker  cannot  be  reduced  to  simple 
expressions  in  the  variables  and  in  the  parameters  as  is  the  case  for  functions 
of  the  hypergeometric  type;  in  fact  these  coefficients  involve  products  of  the 
functions  themselves.  Therefore,  although  Whittaker 's  formulas  for  the  Mathieu 
functions  have  the  form  of  recurrence  relations  they  do  not  have  the  properties 
of  recurrence  relations  for  fimctions  of  the  hypergeometric  type. 

The  procedure  used  in  the  present  paper  is  similar  to  that  used  by 
Meixner  and  Schafke'-  -■  for  Mathieu  functions.  The  different  integral  rela- 
tions are  obtained  by  choosing  different  kernels  and  taking  various  solu- 
tions of  the  spheroidal  equations.  Among  the  relations  thus  obtained  some 
3d.eld  just  the  coefficients  of  the  'recurrence  relations'  for  spheroidal 
functions.  As  in  the  case  of  Mathieu  functions,  these  coefficients  cannot 
be  reduced  to  expressions  which  do  not  involve  products  of  the  functions 
themselves,  Marx^  -'  obtained  some  of  the  integral  illations  derived  here. 
However,  our  procedure  enables  us  to  avoid  the  use  of  expansion  theorems 
for  spheroidal  functions,  and  our  results  are  more  general. 

We  have  also  investigated  the  limiting  case,  when  y  =  0*  which 
leads  to  a  number  of  new  integral  relations  involving  products  of  associated 
Legendre  functions. 
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2»     The  wave  equation  in  spheroidal  coordinates 

The  three-dimensional  time-harmonic  wave  equation 


(2.1) 


Au  +  k  u  =  0, 


where  k  is  a  constant,  can  be  separated  in  prolate  spheroidal  coordinates. 
These  are  defined  by 


(2.2) 


X  -  c  /(4^-l)(l-?}^)   cos  0, 
y  -  c  7(4^-1  )(1-;|  2)   sin  0, 


z  «  c  4 


1 


If  we  take  a  solution  of   (2.1)  of  the  form 


(2.3) 


u(x,y,z)  =  U(C)  V(7[)  W(0), 


then  the  wave  equation  reduces  to  three  ordinary  differential  equations 


(2.Ua) 


(2.i4b) 


[<^-^>  \],  ♦  f-  i 


+  Y^(l-4^) 


U  »  0, 


[(i-;,2)v^]    .    [^-_jl-.y2(^_^2)   ,,q^ 


(2.l4C) 


^^^*^^V'0, 


where 


(2.5) 


2    2  2 
Y  -  c  k'^ 


and  [I  and  X  are  the  separation  parameters, 


-  u  - 


The  differential  equations  (2.Ua)  and  (2.Ub)  are  identical^  they  are 
called  the  differential  equations  of  spheroidal  functions.  They  have  been  dealt 
with  in  great  detail  in  the  treatise  by  Meixner  and  Sch'ifke  mentioned  above. 

We  have  the  following  theorem: 

Theorem  I.   If  |   is  a  path  in  the  complex  71   -plane  and  D  is  a 
domain  in  the  5-plane  and  if  ^2^^^   ^"*^  u^(0)  ^^®  solutions  of 
(2.Ub)  and  (2.Uc)  respectively,  and  furthermore  if 

(A)  u(4,7|,0)  -  G(4,>|)  u^CjZf) 

is  a  regular  analytic  solution  of  (2,1)  in  4  and  ?|  in  their  respec- 
tive domains,  then 

(B)  u^(C)  -  r  G(4,7])  UgC^)  d?| 

r 

is  a  solution  of   (2.Ua)   in  D,   provided  that 

(C)  (1-^  2)   f°%p    U3(7,)  -  GU,^  )  f^ 

takes  the  same  values  at  both  ends  of  the  path  of  integration  P 
and  is  such  that  the  indefinite  integral  of  (B)  is  uniformly  con- 
vergent for  all  4  in  D.* 

Many  known  integral  relations  between  spheroidal  functions  follow  in  a 
straightforward  way  from  this  theorem  if  one  simply  chooses  for  u(4,'?7,0)  ele- 
mentary solutions  of  the  wave  equation  with  angular  dependence  e  ^.  Of  course 
one  could  use  the  separated  solution  (2.3)  and  take  as  a  kernel  G(4,  ??)  in  (B) 

(2.6)  G(4,>/)  -  u(C,?^,0)  e"^^^-  U(4)  ViTJ), 

A  simple  proof  is  given  in  section  1.133  of  ref .  [U]. 


-  t^  - 


where  U(4)  and  V(>^)  are  solutions  of  (2.Ua)  and  (2<lib)j  but  then  the  integral 
relation  turns  out  to  be  trivial.  This  is,  however,  not  the  case  if  we  take 
Lu  instead  of  u,  where  L  is  an  operator  which  comrautes  with  A  and  is  such 
that  it  produces  an  angular  dependence  e      ^,     Then  Lu  is  also  a  solution  of 
the  wave  equation,  and  by  taking  e"  ^  "^  Lu  as  a  kernel  in  (B)  one  gets  non- 
trivial  integral  relations  involving  products  of  spheroidal  f\inctions. 

Certain  combinations  of  the  components  of  the  linear  and  angular 
momentum  operators  have  the  property  assigned  to  the  operator  L»  These  are 

4  ax  *  ^  ay  ' 
■^z  "  a?  ' 


(2.7)  A      ^  A      ;» 

M^-M  ±iM  .(y|.-z|-)±i(z|--x|-), 

±   X    y  ^"^  az   ay    ^  ax   az  ' 

M  «  (x  —  -  V  — )  . 

z   ^^  ay  ^  ax-*  • 

Powers  and  products  of  these  operators  possess,  of  course,  the  same  property 
as  the  operators  themselves. 

We  also  have  the  following  lemma. 

Lemma  I.   Let  sJJ^-^^C^)  and  sJj/'^^C)   (j-1, 2, 3, U)  be  solutions  of 
the  spheroidal  equation  (2.i4a)  f   then  for  X  =  X^(y  )  and  X  »  xj^i  (y  ) 
respectively  there  exist  relations  of  the  type 


(2.8)     a(4)  S^}^\^)  *  b(4)  S^}^^\g)   -  S^:'.(J^4), 


* 


For  the  definition  of  the  functions  S^^'^^Cc)  see  ref .  {I],  sec,  3.65,  Any 


V 

two  of  these  functions  are  linearly  independent. 

For  v-v  =  integer,  p.-^  =  integer,  these  relations  can  be  considered  as  re- 
currence relations  between  spheroidal  functions;  this  case  has  been  treated 
in  detail  by  Marx. 


-  6  - 

(2.9)  c(4)  s;;/J^4)  +  dU)  S^/«^^4)  -  S^^J^4), 

valid  for  J  ■  1,2,3,U.     The  coefficients  a(4), d(4)  are  given  by 

(2.10)  a(4)  -    [^',v';   y.,   v]*,         b(c)  -    [^l'v';   \i,   v] 

(2.11)  c(4)  -  [ti,  v;  ii',  v']^   d(4)  -  [ti,  v;  ^'v']  . 
The  expressions  In  brackets  in  (2.10)  and  (2.11)  stand  for 

r  .  .    1*   Sj'(^^0s^^^^'(4)-Sj^^^4)si;(^^'(4) 
(2,12a)   [vi  ,v   ;ix,  vj  -  jjYj * 

W  '      (K) 


r ,  .    1       s;;,'(J\4)s^^^^(4)-s5,'(^^4)s;;(^^c) 

(2.12b)   [n  ,v  ;ii,vj   -  ^^j-j^^ 

W  '   (?) 

and 

r      .  n*    s;;('5)(,)sR;(ic)\,).3.(k)(^)3j;;(j)\^) 

(2.13a)   [\i,   V|n  ,v  j  -  ^-^ 


Vv  (^^ 


.      ,  n       s;;^j)(os!;.^^\4)-3^,^"^?)s!;/^^o 

(2.13b)   [ti,V}n  ,v  J   -  -^ ^^ ' 

.    ,    (J.l^),  , 
These  expressions  are  independent  of  |j.,  vj  jj.  ,  v  .  W     ($)  is  an 

abbreviation  for  the  Wronskian 

St.(j)(^)S|X(k)\^^.Sjt(k)(^^g^.(j)\^)    . 
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The  proof  is  obvious.  Furthermore  it  follows  easily  that  no  other 
relations  of  the  form  (2.8)  and  (2.9)  exist. 

The  operators  (2.7)  applied  to  the  solution  (2.3)  of  the  wave  equa- 
tion lead  to  the  following  kernels: 


(2.1U) 

(2.15) 
(2.16) 


We  restrict  ourselves  in  the  following  to  spheroidal  functions  with  j  =  3,U; 
this  implies  no  loss  of  generality  because  the  other  cases,  (j  ■  1,2)  can  be 
obtained  by  employing  the  relations 

(2.17)        s^J^^'^h^)  -  s^^^ho  *  i  s^;(2)(u  . 

V  V  V 

The  operators  (2.7)  may  be  written  in  spheroidal  coordinates  9S 


(2.18) 


(2.19) 


z   c 


2  2 

4  -  1  ^  3  +  1-^^  _  a 


\ 


■»•   c 


'«  as  -  1  5^  '  ^ 


2   2 


ie 


±10 


7(4^-1)  (1-ri  2) 


a 
a^ 
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(2.20) 


M. 


^^^  ink-^k''''--?=^ 


^  '^  '      ViC'Dii-^ 


For  these  three  operators  respectively  we  take  as  functions  Up(/|  ) 


\/(?^i)a^ 


(2.21)      ^2^^^'rl2p^i(^^  * 


(2.22) 


(2.23) 


(nn) 

Up(77)  -^     (-?/), 
'      v+2p+l  ' 


(uti) 


(nil) 


The  functions  J^{?)  )j  511     (>/)>  etc.,  as  well  as  the  functions  s|J^^  , 

are  solutions  of  the  spheroidal  equations  with  X  »  '^^(y  )>'*•"  ^^^2d+1  ^^  ^» 
etc.  Here  one  should  remark  that  p  is  an  arbitrary  nuraber^  the  notation 
V  ■  v+2p+l  and  v  »  v+2p  is  convenient  because  in  the  special  case  (A  )  (see 
below)  one  has  to  choose  p  as  an  integer. 

The  path  of  integration  \  in  (b)  is  chosen  in  such  a  way  that  the 

condition  (C)  of  Theorem  I  is  satisfied.   This  choice  may,  of  co\arse,  also 


impose  conditions  on  the  admissible  functions  J]^   (■??)•   We  then  arrive 
at  the  integral  relations  of  spheroidal  functions 
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(2.25)  C  4'^rs^X-^^^^'  ^ 


■'n  r-^  '         V+2D+1  '       ^ 


T 


;ja)(,)|      r/(^^-l)(l-^^)^q'(^) 


V^      ^'1  '     ^        s/a'-Dd-n') 


(titl) 
E  (^)d?7 

v+2p+l     '         ' 


V        I 


I        ,  I 

where  |x  ■  till,  v  «  v+2p+lj  and 


%'  3         v+2p         ^^^ 


(2.26) 


±± 


ij"''(«)j   ^SpnpS^  £:i^(^ 


v+2p+l    ' 


where  r  -  ixt  1,  v  -  v+2p,  and  where  a^,  p^,  and  6^,  are  constant  quantities 
determined  from  the  asymptotic  behavior  of  the  kernels  G  ,  G  and  G. 
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[eqs.    (2,lU),    (2.15)   and    (2.l6)].     One  obtains  thus 


t 

(2.28)    ■  f'\    J      A'Tf     T^{7[)Y.  (^')d7, 

V  V    I         v+2p+l    I         ^ 


r       ft o  '  \  (^^±1) 


(2.29)  _ 


The  e^*'  ,  (J»l,2,3,li),  (k"l,2,3),  are  phase  factors  depending  on  v  ,  v  of  the 
spheroidal  functions  S^^''^(c).  For  instance  if  we  set  (J-3,i4),  we  find 

(2.30)      e^^'*^^  -  e^i(P*l)«      (k.1,2,3) 

where  the  plus  sign  is  taken  when  j«3  and  the  minus  sign  when  j-U.  The  phase 
factors  e^  '   can  be  determined  from  the  asymptotic  behavior  of  ^  *      (4) 
or  directly  from  the  relations  (2.17). 

Application  of  Lemma  I  leads  then  immediately  to  the  following  integral 
relations  involving  products  of  spheroidal  functions: 


-  u  - 


(a.,3)  1/     ^SpS  ,  j:.,^,  ^^\^  ,  ,^   .  Pvap)[.,v„.l,v,.p4 


p  4-^  \        v+2p' 


PL     '"^         '      v^(?^ 


prep^Lil,   v+2p+lj  ^i,vj*  , 


,,„         -  r   ^^^^^^Zy  iT^  a^  -  a  eO)[.v..n,.a.]. 


n  4  -^  (      V  I        VH-2P+1 


(2.36)      T i /  W-^nyf^ ,  c:"' (^)  t     ^y      v^^-  - '  -^^""^ 
rL    ^"^        "  ^   /(4^-i)(i- 


•^') 


^^:<'(Y^..p^'^'[ 


* 


.  6  e^J^ 


[|itl,   v+2p,-   ti,v|. 


where  a  -  aJJ,   P  -  pj),   5  -  sJJ  . 
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3»  Special  kernels  , 

The  resTilts  of  Section  2  are  general  in  so  far  as  they  are  valid  for 
any  function  ^^{^  )  which  is  a  solution  of  the  spheroidal  equation  and  for 
a  path  P  satisfying  condition  (C)  of  Theoran  I.  Now,  we  shall  consider 
specific  solutions  of  the  spheroidal  equation  by  restricting  the  path  of  in- 
tegration and  the  parameters  m.,  v  and  y. 


Case  A.  Let  (v-|x)  be  an  integer  p.  The  spheroidal  functions  y^(Y\  )  then 

become  ^^(?7)  -  psj;(7l  JY  ),  iZ         ("l^  '^Vap+l  ^^  '^  ^*   ®*'^'  ^  appropriate 
path  of  integration  P  is  then  the  real  interval  [-l,l].  We  should  remark, 
however,  that  ps^(7J ;y  )  remains  finite  at  ± 1  if 

Re  |x  <  0      for       V  +  ji  -  0,  1,  2,'»« 

or 

Re  n  <  0      for       V  -  ti  -  -1,  -2,  -3,*»»  . 

If  Re  iJ.  >  0,  then  p8^(?7  ff  )   will  in  general  be  infinite  for  77  -  +  1, 
except  when  \i  and  v  are  both  integers,  in  which  case  it  remains  finite  at 
>?  «  ± 1,  In  the  latter  cas^  the  conditions  of  Theorem  I  are  fulfilled,  as 
one  sees  from  the  fact  that  the  functions  ps^(?|  JY  )  are  then  even  for  v-n  even 
or  odd  for  v-p,  odd,  and  behave  as  (1-7))°''   in  the  neighborhood  of  ??  «  tl. 
But  there  is  no  loss  of  generality  from  the  restriction  Re  ix  <  0,  because  for 
Re  M-  >  0  there  exists  a  certain  combination  of  ps^(  W  ,Y  )  and  qs^(  hfr  )   ^ich 
is  finite  for  v-p.  «  0,  1,  2, ...or  v+n  ■  -1,  -2,  -3...,  and  then  it  can  be  ex- 
pressed by  ps"^(y1;Y^)  with  Re  ti  >  0,  so  that  this  case  is  reduced  to  the  former 
one. 
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It  is  fxirther  necessary  to  distingtdsh  between  the  cases  (2»2U)j 
(2.25)  and  (2.26)  and  one  must  check  in  each  particular  case  whether  the 
condition  (C)  of  Theorem  I  is  satisfied. 

Case  B.  Let  J^{y\)   «  ps'^iYf  fr  )»     For  real  y  >  0,  we  take  for  P 

V 

the  interval  [l, ..,oo]|  otherwise,  we  have  as  an  appropriate  path  P 

1, ...00  arg  -•   for   Re  y  >  0 
and 

1, ...  00  arg  Y   for   Re  y  <  0  , 

It 
where  y  is  the  conjugate  complex  of  y. 

Case  C.  For  ^  ( ??  )  arbitrai^r,  we  choose  P  to  be  a  Jordan-Pochharamer 

V 

path  (double  loop  around  any  pair  of  singular  points  of  the  spheroidal  equation). 
Case  D.  If  we  take  ^^^(>|  )  -  sJJ^^^?])  or  SJJ^^^)]),  then  the  argument 

V 

of  ?1    for  Ty  ->oo  has  to  be  chosen  in  such  a  way  that 

|arg  (-iyOl  >  |   or   |arg  (iY4)|  <  ^ 

respectively.  The  path  of  integration  is  then  a  single  loop  starting  from 
ioo  e   ,  passing  around  the  point  +1  in  the  negative  direction,  and  return- 
ing to  infinity  in  the  direction  oo.e^^"^"^^/^.  For  S^^^\rj)   and  for 

SJ^^^A]  )  the  path  P     starts  at  oo-e^^^""^^/^  and  after  passing  around 

in 

2  "" 

+1  in  the  negative  direction  returns  to  infinity  along  a  line  oo.e      . 

(See  Figs.  1  and  2.) 
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Figure  1 


Figure  2 


Case  E«  If  we  choose  J^^iT))   -  sJJ^^^C)^  )  then  P  is  a  loop  similar 

V 

to  that  in  case  D,  but  now  the  path  goes  round  both  points  -1  and  +1  (see  Fig.  3). 


Figure  3 


-  15  - 


f 

Case  A  «  An  important  instance  of  case  A  occurs  when  we  restrict  p, 

and  V  to  integral  values,  v  ■  n,  p.  ■  m  and  n  >  m  >  0.  In  this  case  the 
spheroidal  functions  are  either 

(mtl)        (mil)    2 


^n«p'*^'"''"*2p'*?"' 


or 


(rail)         -(mia)     „        (mil)        -(mil) 


Thus  we  can  expand  the  spheroidal  ftmctions  ps   ()9jY  )»   etc.,   in  associated 
Legendre  functions  of  the  form 


(3.1)  KC/'^')-    f^    <,2r(v')  CarO?  >. 


»-2)  P'nVl'1 '^'  -  J^    "nVLarfv")  C2p»>2r(>?) 


The  constants  a,  p,  6  are  given  by  the  expressions  (2.27)- (2.29), 
which  can  be  easily  evaluated  from  the  properties  of  associated  Legendre  func- 
tions. For  instance,  we  find  the  following  expression  for  the  constant  o  -  a"  (r^)- 
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a  -  a"*  (y2)  -  f^-y 


c(2n+Ur+i) 


oo 


r»-oo 


m     2)  -m         ,  2x  n^2r-m^l 
n,2r^^  -^  n+2p+l,2r-2p^'  '   2n+Ur+3 


(3.3) 


i-^      n,2r''^  ''  n+2p+l,2r-2p-2^''  ^  2n+Ur+l 


The  other  constants  p,  5  can  be  found  in  a  similar  way  by  substituting  for 

u     2     '^^'^'^^ 
the  spheroidal  functions  VK,^Y)  if   )»  PSv+2d+1»  ®**^»»  "^^^  corresponding  ex- 
pansions in  associated  Legendre  functions  in  (2.28)  and  (2.29)  and  then  per- 
forming terniwise  integration. 


k»    Limiting  cases 

There  are  two  limiting  cases  for  which  the  spheroidal  functions  re- 
duce to  simpler  functions.  If  we  put  Y  "  0»  then  the  differential  equation 
(2.Ua)  goes  over  into  the  differential  equation  of  the  associated  Legendre 
fvinctions.  We  therefore  expect  that  our  integral  relations  reduce  in  this 
limiting  case  to  integral  relations  between  products  of  associated  Legendre 
functions  of  the  first  and  second  kind. 

However,  if  we  introduce  y4  "  z  as  a  new  variable  and  then  set  y  "  0, 
the  integrands  in  (2.2U),  (2.25)  and  (2.26)  split  into  two  factors,  one  de- 
pending on  z  and  the  other  on  yj   .     Therefore  this  limiting  case  of  our  in- 
tegral relations  gives  a  rather  trivial  results 

We  now  shall  proceed  to  work  out  these  limiting  cases  in  detail. 

Case  I.  Y  "  0.  It  is  easier  to  treat  the  limiting  case  in  a  direct 
way  instead  of  working  out  the  limiting  process  from  our  general  results 
arrived  at  in  the  previous  sections  .  Since 


This  alternative  method  is  given  in  Appendix  A, 
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Ps"(>j^|Y^)  -►  P^()^)     for     Y-0 


while  S^^^*^^'Ufr)   reduces  to  a.  linear  combination  of  ?^^U)   and  Q^^C^), 
we  can  at  once  write  down  the  following  formulas 


m  ^tM  0  t  m 

-1    ^  "   1  ' 


-1  ^  "  ( 


and 

1 


-1  ^  "U 


where  the  coefficients  a  ,  a_ j  a  ,  a  are  to  be  determined.  From  the  behavior 
for  large  i   of  the  integrals  and  of  the  functions  on  the  left  side  of  (U.l) 
and  (U.2)  one  immediately  concludes  that  a.  ■  0  sind  a  -  0  and  furthennore  that 


(U.3)       a^  -  0    for     P  <  0,     a    '  0        for    p  >  0. 


Similar  expressions  can  be  written  for  P°(4)  and  (^"^(4)   (see  remark  following 
Eq,    (U.18)). 
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In  the  other  cases  one  finds,  expanding  in  descending  powers  of  4  and  com- 
paring the  coefficients  of  the  leading  terms, 

(U.U)       a^  -  2(-l)"*^  (p  >  0),    a^  -  2(-ir  (p  <  0)  , 
Solving  for  the  two  integrals  occurring  in  (U.l)  and  (U»2),  we  get 


i-1 


-in  '  -m  ' 


(i"*'  "l  «n.2p*l(«)  ■"„  («>  -  "2  '■„.2p»l(C)  «:  (4) 


"'■""J        X&    ^n'('l)C2p.l(1)^n 


--1      «-1 
In  a  similar  way  we  get  from  (2.25)  in  the  limiting  case 


f    -(mtl)  ml    -(mil) 

(»^-7)  Pi"  V2p.l(^^  *  ^1     Pn.2p.l(«) 


/^         /f^2  ,s/,     u2n  ■  •  -(mH) 

J        ^^y^'^'-?    '       «  C  <«'  ^^C?  '  fn*2pn(7  >  "V 

-1       ^  "'( 

-/        ^%''^T1    ^       n    C(^)  Pn   (^)  ^n.2p.l(7^  '1 
-1  ^  "^' 

/(«^-l)(l-^2) 
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mt    -(mtl)  mt    -(mtl) 


■  /  "^"^^'y^y'^ '  '""'^'^  "^^"^  ^  '-"^'^^"^ ^  '"^ 


-1  ^-^ 


r^       /rr2  T\n    h2;  „  „'  -(mil) 

-1  ^  "M 


while   (2.26)  reduces  to 

+    -(mtl)  rat    -(in±l) 

(U.9)     6^;-  Q„,2p.l(C)  *  5^     Pn.2p.l(^) 


,,J      /^^'-jM>>?^)'    v)  q;-'(,)  p;(;^  )  i^^\yi)  d/; 


'-1        ^-^ 


_,j   /(^i)(i->^^)'  ^ ,-.(^) p.;^) p'j^r^o/) ^>? 


'-1      ^-^ 


-(mil) 


-^  {     ^  ,   '^     ..  «;"(*'  ^(1 '  Cap'   (>! '  M 
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and 

mt    -(mtl)  _ni±    -(^l±l),.^ 


/     /iTTTuTT^        J       m        -(rati) 
t,j     A^  -i)a;t  )  >,  P-  (o  p°(>, )  p„,,p    (/,)  d/i 


c-^ 


Following  the  same  argxunents  as  for  the  a     and  p.    (j»l,2),  we  find 
m±  ^+  mt  m-j  ^  *" 

p       »  0  and  pp     -  0  and  5,-0  and  5.     ■  0.     Furthermore,   one  has 


m±  m± 


(h.ll)  3,-0       for       p  <  0,  P2     "  °       ^°^       p  >  0 


and 


ml  m± 


(U.12)  6^-0      foi-       p  <  0,  62-0      for       p  >  0 
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In  all  other  cases  one  finds,   expanding  in .descending  powers  of  4  and  comparing 
the  coefficients  of  the  leading  terms. 


,m+  n-»it>+2  /  ,  \ra+l  x  n-Ei-1 


(U.13a)  Pi     '  ■  2ivia~*^'^^       '         P 2     ■  " 


2n+l 


^m-  n-m+2  /  ,  \m-l  7^-  "  n+m-1 


(u.i3b)        p^  -    2^^^(-i)    ,     P2  ■    2H7r 


/,    -i„\  ^m+       ,  n+m+2  /  ,  \m+l  r  .  n-m-1 

(UolUa)  5^     -  i  2H5j-(-l)       ,         62  -  a  ^^^ 


f\    iiv.\  R^~  _     *   n-m+2/  T  sin-l        —  ,   n+m-1 

(U.lUb)  5^     -  -i  2H^(-1)       ,       62  --  i  5^^  . 


Solving  for  the  two  integrals  occurring  in   (U.?)  and   (h.8),   we  get 


mt    -(nttl)  ml    -(mtl) 


1       r^    X~«7      „        -(""ii) 
5=  /^    ^  «0l>^„.ap.i(1>^'). 


^C^-1 
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mt    -(mil)  '  mt    -(mtl) 


*/. 


^n+2p+l 


(•|)d1. 


and  from  (U.9)  and  (U.IO)  one  obtains 


mt    -(mtL)  mt    -(mil) 


-(mtl) 


nOi'Vap  'vn 


and 


mt    -(mil)  J  m      -(mtl)  » 

^1    «n»2p     («'  "-n  (4)  -  6^    P„,2p     (U  Q:  (4) 


(U.18) 


,^J^(m^^\^,. 


m 


4 


/(4-l)(l-/7^) 


^n^V 


-(mtl) 
^n+2p 


T 


(^)  d^ 


Finally  we  wish  to  remark  that  if  we  replace  m  by  -m  and  -(mtl)  by 
(m+1)  in  the  associated  Legendre  functions  one  gets  similar  integral  express- 
ions which  involve  associated  Legendre  functions  and  whose  coefficients 

«mt  _«mt   «m+  _»m+ 
P-i  *  Fo  f  ^1     »  ^2         ^^®  obtained  from  those  of  (U.13a),  (U.lUb)  by  replac- 
ing m  by  -m. 
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Case  II.  Let  us  replace  4  by  z/y  and  then  put  y  ■  0»  The  results 

m       O 

of  Sections  2  and  3  are  easily  applied  to  this  case  by  replacing  ps  {Yl  if  ) 
by  f^in)   and  sJJ^^^(4)  by  Y^^^(z).  The  functions  ^l^\z)   are  defined  by 

(J) 

where  for  J  ■  1,2,3,'4  respectively,  Z  ^^lAz)   is  the  spherical  Bessel  function, 

the  Neumann  function,  and  the  Hankel  function  of  the  first  and  the  second 

m    ^ 
kind.  The  remaining  task  is  the  evaluation  of  the  constants  a   ,8   , 

mH  ^'P       ^*P 

and  5  for  Y  =  0,     One  finds 

n,p 


m      ^      2       n-m*l 
°n,o  "  2n+l     2n+l       * 


m  2        n+m 

a,         ~ 


n,-l     2n-l     2n+l   ' 


„m+l  2  1  p,m+l     _  2  1 

Pn,o  '  2n+l     2n+5'  '  '^n,-l       "  2n+l     2n-l  ' 


m-1  2         (n-nn-l)(n-m+2)  qIti-I     ^        2        (n-«-m)(n-<-m-l)  . 

Pn,o  "  '  2n+l  2n+5  '  ^^n-l  2n+T  2n-l  ' 


g^m+1  n-gn-1 

n,o  ■     C2n+inn+m+i;  * 


m-1  _     (n-t-2n-2)(n"ra-»-l) 
n,  o  2n+l 
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All  the  other  coefficients  vanish.  Finally,  the  expressions  (2.12a)  and 
(2,1b)  reduce  for  j  »  3,U  to 

(3)     (U)'      (U)'     (3)'     n[U,3] 
Vl^^^^n   (^)-Vl(^)^n   (^)--T-* 

(3)     (U)'      (U)     (3)'      (n*l)[U,3] 
Vl(^)^n   (^)-Vl(^)^n   (^) 1 


and 


(3)    (U)      (U)    (3)      [U,3] 


(3)    (3)      Ox)  (3)       [U,3] 


where  [3,U]  is  the  Wronskian  T^^z)  T^^^  (z)  -  I^^\z)  T^^^  (z).  These 
last  results  are  rather  trivial. 
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Appendix  I.  Integral  relations  involving  products  of  associated  Legendre 
functions 

In  Section  U  we  derived  integral  relations  between  associated 
Legendre  f\mctlons»   Here  we  shall  give  an  alternative  method  of  deriving 
such  integral  relations  by  considering  the  general  formulas (2. 3li),  (2.36) 
for  Y  small  and  then  taking  the  limit  as  y  -♦■  0» 

In  general  we  can  write  (2.31)  and  (2.32)  in  the  form 


(1.1) 


a*   Y^  e-^^P^l)" 
n,p 


■t(3,U)    ni(U,3)      t(U,3)    in(3,U) 
Sn.2p.l(^^  \  (^>  -  \.2p.l(^^  ^n  ^^^ 


•  ip(^)  t^'3]  , 


(1.2) 


^t    2  gll(p+l)n 
n,p 


■t(3,U)    m(U,3) 

\.2p*l(«^  ^n      (^)  -  Sn.2p.l(^)  ^n  ^^^ 


t(U,3)    ra(3,U) 
^2p^ 


\^U)    [U,3]  , 


,        ,  t(3,U) 

v^ere  p.  -  m,  v  -  n  and  m-  ■  t,  v  -  n+2p+l  or  n  +  2p.  We  can  replace  S     ($) 

t(l) 
S    ($)  using  the  relations 


(1.3) 


t(3) 


t(l) 


-ivn 


t(l) 


cos(vTi)S^  U)   -  -1  S_^_^(4)  +  e"^^"  S^   (S), 


We  denote  the  Wronskian  of  the  spheroidal  functions  by  [3,U]. 
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t(U)      t(i) 


+ivn 


t(l) 


(I.U)      cos(vii)S^  U)   -  i  S_^.^(4)  ♦  e  """  S^   (4) 


We  introduce  the  abbreviation 


(1.5) 


■t(3)      m(U)      t(U)      ra(3) 
Sn.2p.l(^^  ^n   ^^)  -  2n.2p*l(^)  ^n   ^^^ 


t(l) 
or,  in  terms  of  S    (c). 


(1.6)      A  -  2i(-l)" 


"td)       md)      t(l)      m(l) 


Writing 


(1.7) 


*^1^       t     t    2 


t(i) 


-itn 


(1.8)         s.^.i(?^Y)  -  -t — r-7-::r  2, 


e"""cos(vn)n  (v-t+l)      t 


.  QsUfr^) 


r\ir)\ir'')\  ir)r  (v+t+i) 


and  substituting  in  (1.6),  we  find 


(1.9)     A  -  2i(-l)"' 


r  (n-t+l)V°(Y)V*^  pn-?(Y^)     -t       p    m    2 
n  (n+2p+l-t+l) 


r  (n+2p+l+t+l)   t        „        -t       o    „    o 
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where 


(I.IO)      kJ(y)  -  cos(vn)  P  (v-t+1)  vJ(y) 


and 


(i.n)     Y  v^(y)  v.^_^(y)  a*(y^)  a^  (y^)  -  1 


For  small  values  of  y  the  function  V  (y)  is  given  by  the  expression 


(1.12)      vJ(y)  -  I  (iJ       1^  (l  +  0  (Y^))  . 

r   (v-^  |) 

Substituting  this  expression  into   (1.9)  we  obtain  for  small  values  of  y 


(1.13)  A  -  2i(-l) 


m 


r  (n-t+l)  r  A  -n)  r  (n+2p+  |)      ^  -2p-2 


h   (r  (n+2p+l-t+l)  r  (n+  |)  r  (-n-2p-  |)) 


(f) 


^V2p.i(^^^'^  K^^f'''^ 


r  (n+2p+l+t+l)  r  (n+  |)  r  (-n-2p-  |)    ^  2p 

— 7 \ —  %^ 

U(r  (nn+l)  r  (|  -n)  P  (n+2p+  |)j 


^Wl^^^^')  ^^n^^^^'^ 


See  ref.[U],  page  297. 
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and 


(I.lU)  B  -  (2i)(-l) 


m 


r  (n-m+l)  r  (|  -n)  P  (n+2p+  |)        -2p-2 


'— _£ : _i (J) 

U  (r  (n*2p+l-t+l)  r  (n*  |)  r  (-n-2p-  |)) 


«C2p.i(^^y'>  Psf  (4;y') 


r  (n+2p-H+t+l)  r  (n+  h    r  (-n-2p-  i) 


U  (r  (nH+1)  r  (|  -n)  r  (n+2p+  |)j 


* 


2p 


pC2p*i(^'^'^  Q3^'(4^y2) 


where  B  stands  for 


t(3)      ni(U) 
Sn.2p.l(^^  ^     (^^ 


t(U)      m(3) 
Sn.2p.l(^>  ^n    («^ 


Therefore   (I.l)  and  (1.2)  become 


(1.15) 


ip  re'^'^"'^  A     .     .IJO    [U,3]   , 


(1.16) 


a*       Y6^(p+l)n         .       t      (^    n      -J 


vd-th  A  and  B  given  by  the  expressions  (1.13)  and  (I.lU)  respectively. 
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Now  we  are  in  a  position  to  take  the  limit  y  ■  0«  The  spheroidal 
function  Qs  (C;y  )  and  Ps  (CjY  )  go  over  to  the  associated  Legendre  functions 
Q^  U)   and  P"(4).  Hence  (3.1U)  and  (3.15)  reduce  to 


(1.17)    lim 

Y  ->  0 


CnYVi^P*'^"(2i)(-l)"' 


,     r  (n-«-^l)  r  (|-  -n)  r  (n+2p+  |) 
U  (r  (n+2p+l-m+l)  r  (n+  |)  P  (-n-2p-  |)) 


-2p-2 


•  C2pn(«>  ^>^ 


r  (n*2p*l*nn)  r  (n*  |)  f  (-n-2i>-  J)    ^  2p 

-. «_  (^j 

U  (  r  (n+m+l)  r  (|  -n)  P  (n+2p+  |)j 


-m 


Dl/ 


(4)  Qr(4) 


n+2p+l^^'  n 


-1  ^  ~  I 
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and 


(1.18)       lijn 


^m      YVi(P*l)«  2i(-l)" 
n,p 


r  (n-m+1)    r  (|  -n)"   P  (n*2p+  |)  -2p-2 

L  ('r  (n+2p+l-m+l)    r  (n+  |)    P  (-n-2p-  |)) 

I 
-ra  +m 


r  (n+2p+l+m+l)    P  (n*  |)    P  (-n-2p-  |)       ^  ^P 
U  fp  (n+it+1)    P  (|  -n)    P  (n+2p+  |)) 


(f) 


••m  n 


fifr  /  piV  ^  <i '  Cpn  ^1 


r-i    -Li  r-1 


It  should  be  remarked  that  a"   is  a  function  of  Y^,  as  can  be  seen  from 

n,  p 

Eq.  (3.7)  or  (3.11).  It  is  evident  that  the  coefficients  aj,2r^^  ^  ^*^ 


a 


„  ,  o  .  (y^)  are  proportional  to  r^^   and  Y  '  >  hence  the  first  terms 
n+2p+l,2r-2p'''  '  ^     ^ 

of  the  left-hand  sides  of  (1.17)  and  (1.18)  are  independent  of  Y.  The  same 

argument  applies  to  the  second  terms  of  these  expressions.  It  suffices  to 

show  that  (1.17)  and  (1.18)  agree  with  the  results  derived  in  Section  U  for 

p  «  0  and  p  ■  -1.  In  both  these  cases  a    is  independent  of  Y»  t-h^tt  is, 

n,  p 

<^     {i)   ■  '^t   and  is  given  by  the  integral 
n,p 
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-1 


We  get 


m    in       m     n  (n-m+l)  P  (i  -n)  r  (n+  |)      -m 

(1.20)  a^^^  e^\2±)i-ir    ^ L~  U  Q^^l^^^  ^n^^^ 

r  (n+2p+l-m+l)  r  (n+  |)  R-n-  |) 


-1   ^   ' 


(1.21)  a»    e^"(21)(-l)"  ^ U  P„,^(4)  Q»(4) 

r  (nwl)  r  (t  -n)  Rn*  i) 


f  /  75^  Ol^ClO''"? 


To  prove  that  (1.20)  and  (1.21)  are  the  correct  integral  relations  be- 
tween products  of  associated  Legendre  functions  one  must  show  that  for  large  4 
they  reduce  to  identities  in  4.  For  this  purpose  we  have  only  to  consider  the 
leading  terms  of  Q  (4),  r  (4),  etc.,  and  compare  the  coefficients  of  both  sides 
of  (1.20)  and  (1.21)  in  4.  Then  the  statement  is  easily  verified*. 

In  the  second  limiting  case,  Y  ■  0  —  ■  4,  the  results  are  rather  trivial, 
since  the  constants  a    vanish  for  all  p  /  0  and  -1  (see  Sect.  3,  Eq«  (3.7)). 
When  we  put  p  =  0  and  p  ■  -1  and  remember  the  limiting  expressions  for  S    (4,Y  ) 
and  ps^^iij  ir  ),   we  find  from  (3.12) 

We  remark  that  for  p  =  0,  only  the  first  terms  of  (1.17)  and  (I. 18)  give  a 
non-vanishing  result,  and  for  p  ■  -1,  only  the  second  terms  remadn^as  can 
be  easily  seen. 
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(5)  (J)'     „  (J) 

(1.21)  -Y^,,(z)  -  f^   (z)  .  H  Y^  (^)  , 


(j)       (j)      /_+^^   0) 
(1.23)  -Vl(^)-\  (-)*^^n  (^>' 


(j) 

( 
n 

relations,  we  find,  respectively. 


where  T   (z)  is  defined  in  (U.19).  By  adding  and  subtracting  the  above 


(J)         (J)      o„+T    (J) 
Y  ,(z)  +  Y  ^t(z)  -  ?2li  T   (z) 
n-1^       n+1^      z    n  ^  'j 


(j)  (J)         (J)'       T   (J) 

Vl(^>  -  ^n.l(^)  ■  2  ^n   (^>  4  ^n  ^^^^ 


which  are  the  well-known  recurrence  relations  between  three  contiguous  cylindrical 
functions  ■-  -^  . 
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Appendix  II.  Useful  formulas  involving  associated  Legendre  functions 

We  list  here  the  following  formulas,  which  have  been  employed  in  the 


texts 


(II.l) 


(II.2) 


J 7   ™       T     /n»+l       m+1    \ 

/iW]^  P"(K?)-2H7r  ((n-ni+l)(n-m+2)  P^^3^(/j)-(n+m)(n+m-l)P^_^(;^)) 


(II.3)      /l-^2  Pn  (n>  -  I  [(n+m)(n-m*l)  P^'^C^])  -  P"*^(r|)] 

(II.U)      (l-/]2)  P^;  (/|)  -  ^  [(n+l)(n+m)  P;;.i(/^)  -n(n-m+l)  P^^^C^/)] 


'   n   '     /   A   n  I      n 


(II.6)      /wr  Pn  (^  )  ♦  -^!^  ^(^  ^  ■  ("-^^  ("-'"■'^^  ^n 

/Tip 


ra-1 


2x  ^m, 


From  (II.l)  one  can  write  immediately  the  expression  for  (1-^  )  P  {.T[)i 


(II.7) 


(i-/)2)  p^O]) 


m+2 


m+2 


(2n+l)(2n+5)  V2^V  "  {2n+3)2n-l  ^n  ^V  ^ 


m+2 


2n+ 


^  2^  Vl^V 
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and  by  induction  one  obtains  the  general  case 


q^  (2n-2q-l)i:(2n+2q-Uq+l)(P)       m+p 


q 

^rfiere 

(2n)li  -     (2n+l)  (2n-l)  (2n- 3)   ...  , 
and  where  frcan  (II  .2) 

o  D    n,  SlE          q     (n  +  in)l(n+2p-2q-ra)' 

(II.9)    (i-f)"^  P^(7)  -  IT  (-1)     -; -7; —7- 

^    '  q^  (n-Bi)I(n  +  m-2q)I 


(2n-2q-l)::  s     "»-P         ,^n 

(2n.2q.Uq-l)P^,p.2q('?) 


(2n  +  2p-2q  +  l):: 
From  (II. 8)  and  (II. 9)    one  can  easily  evalroate  the  integrals 

(11.10)  f  f  z:?'  i^c^)  p^"  (?>*'? 

-'-I 

-       0  for  q  /  2p 

22p+2    p(2n+l)(n  +  2p+l);(n+2p-m): 
nl  (n-ra):  (2n+Up  +  3): 

(n.n)  f'f/^p;"(^)p;!;?,'('?>^'; 

"   0 

_22p+2  p(2n)(n  +  2p  +  l):(n+2p-m)I 


for  q  =  2p 


(n-l)l  (n-m):(2n  +  Iip+3): 


for  q  /  2p+l 
for  q  «  2p+l 


With  the  aid  of  (11,8)  and  (II. 9)  one  can  evaluate  integrals  involving  pro- 
ducts of  three  associated  Legendre  functions. 
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